Solute trapping and diffusionless solidification in a binary system 
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Numerous experimental data on the rapid solidification of binary systems exfiibit tlie formation 
of metastable solid phases with the initial (nominal) chemical composition. This fact is explained 
by complete solute trapping leading to diffusionless (chemically partitionless) solidification at 
a finite growth velocity of crystals. Special attention is paid to developing a model of rapid 
solidification which describes a transition from chemically partitioned to diffusionless growth of 
crystals. Analytical treatments lead to the condition for complete solute trapping which directly 
follows from the analysis of the solute diffusion around the solid-liquid interface and atomic 
attachment and detachment at the interface. The resulting equations for the fiux balance at 
the interface take into account two kinetic parameters: diffusion speed Vdi on the interface and 
diffusion speed Vd in bulk phases. The model describes experimental data on nonequilibrium solute 
partitioning in solidification of Si-As alloys [M.J. Aziz et al., J. Cryst. Growth 148, 172 (1995); 
Acta Mater. 48, 4797 (2000)] for the whole range of solidification velocity investigated. 



PACS numbers: 81.10.Aj; 05.70.Fh; 05.70.Ln; 81.30.Fb 



I. INTRODUCTION 

The concept of "solute trapping" has been introduced 
to define the processes of solute redistribution at the in- 
terface which are accompanied by (i) the increasing of 
the chemical potential [1] and (m) the deviation of the 
partition coefficient for solute distribution towards unity 
from its equilibrium value (independently of the sign of 
the chemical potential) 

In experimental investigations of rapid solidification, a 
complete solute trapping leading to diffusionless (chem- 
ically partitionless) solidification was first observed by 
Olsen and Hultgren and Duwez et al. in experiments on 
rapid solidification 0. They showed that rapidly solid- 
ifying alloy systems lead to the originating of supersat- 
urated solid solution with the initial (nominal) chemical 
composition of the alloy. Later on, crystal microstruc- 
turcs with the initial chemical composition were found 
by Biloni and Chalmers in rapidly solidified pre-dendritic 
and dendritic patterns 

Backer and Cahn [l!] have shown that with the finite 
solidification velocity in a Cd-Zn system the coefficient of 
the Cd distribution becomes equal to the unit that char- 
acterizes diffusionless solidification. This fact has been 
confirmed in many binary systems by Miroshnichenko 
0. He investigated dendritic crystal microstructure af- 
ter quenching from the liquid state by splat quench- 
ing and melt spinning methods. The results of Mirosh- 
nichenko's microstructural analysis show that at a cool- 
ing rate greater than some critical value (depending on an 
alloy and experimental method this value is in the range 
10^ — 10^ K/s) a core of main stems of dendrites has initial 
(nominal) chemical composition of the alloy. A critical 
value for undercooling in the transition to purely ther- 
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mally controlled growth with a homogeneous distribution 
of chemical composition in Ni-B solidifying samples pro- 
cessed by an electromagnetic levitation facility has been 
obtained by Eckler et al. Finally, it is necessary 

to note that many eutectic systems undergo chemically 
partitionless solidification with an initial composition [5| 
that can be explained by the transition to diffusionless 
solidification [7]. 

As a consequence, experimental investigations [l|, d, 
0, H, @ show that with increasing driving force of so- 
lidification solute traps are much more pronounced by 
solidifying microstructure. At a finite value of the crit- 
ical governing parameters (undercooling, cooling rate or 
temperature gradient) complete solute trapping occurs. 
Because the finite value of the governing parameter de- 
fines the concrete solidification velocity, complete solute 
trapping and diffusionless solidification begin to proceed 
with a fixed critical growth velocity of crystals. 

The main purpose of the present paper is to describe a 
model for solute trapping and the transition from chem- 
ically partitioned to diffusionless solidification in a bi- 
nary system. Using the local nonequilibrium approach to 
rapid solidification, an analysis of diffusion mass trans- 
port in bulk phases together with conditions of atomic 
attachment and detachment on the solid-liquid interface 
is given. 

The paper is organized as follows. In Sec. |lTl previous 
investigations of solute trapping are shortly reviewed. In 
Sec. mil an analysis of solute diffusion leading to pro- 
nounced solute trapping and complete solute trapping is 
given. The nonequilibrium solute partitioning function 
for atoms on the interface is derived in Sec. IIVI A com- 
parison with previous models and experimental data on 
solidification of binary systems is presented in Sec. [V] 
Finally, in Sec. IVII conclusions of the work are summa- 
rized. 



2 



II. PREVIOUS INVESTIGATIONS 

For the simplest case of an atomic system, let us con- 
sider an isobaric and isothermal binary system (the pres- 
sure P and temperature T are constant) with concentra- 
tion Xa and Xb of atoms A and B, respectively. In this 
article, we denote X as the concentration of the atoms of 
B sort. For a brief overview, we summarize the equilib- 
rium and nonequilibrium solute distribution on the solid- 
liquid interface. 



A. Equilibrium 

In equilibrium, the concentration of atoms X at the 
phase interface is not equal from both sides of the in- 
terface due to the different solubility of atoms in phases. 
During the equilibrium coexistence of phases (gas-solid, 
liquid-solid, gas-liquid) the atoms are distributed along 
the interface in consistency with the diagram of a phase 
state. A difference in atomic concentration in phases 
at the interface can be characterized by the equilibrium 
coefficient /cg of the atomic distribution between phases. 
For equilibrium coexistence of phases (e.g., between crys- 
tal and melt, vapor and crystal, crystal and liquid), the 
coefficient kg can be expressed in the general form 



k,iXL,Xs,T) = ^- ^ 



RT I 



• (1) 



In Eq. ([T]), XI and Xg are the mole fractions of the B 
component in the liquid phase {L) or crystal {S), respec- 
tively, R is the gas constant, and A/j,' is the difference in 
chemical potentials described by 

A^' = A/4 - A/i^,, (2) 

with 

A/is=Mss-MsL. ^^J''A^ tJ'AS- t^AL^ (3) 

where A/i^ and A/^'^ are the driving forces for redistri- 
bution of atoms A and B, respectively, which are defined 
by redistribution potentials /i^ and /i'^ for phases L and 
S. The differences A/i^ and A/i^, Eq. ([3]), define the 
sign of A/i' in Eq. ([2]). For instance, if A/j,' is negative 
(A/i^ > A/i'g), one has k^. < \ - the case of smaller solu- 
bility of atoms B in the phase S in comparison with their 
solubility in the phase L. 

As a general characteristic of phase equilibria in binary 
systems, expression ([1]), together with Eqs. ^ and ([3]), 
is usually considered as a measure of the driving force 
for atomic redistribution at the phase interface. It can 
also be considered as one of the main parameters for the 
construction of the diagrams of a phase state. 



B. Nonequilibrium 

Expressions ([I])-© assume local equilibrium at the in- 
terface, which is a useful approximation for many systems 



transforming at small interface velocities. At a large driv- 
ing force for the interface advancing and with increasing 
of the interface velocity, the local equilibrium is not main- 
tained Therefore, the condition for local interfacial 
equilibrium was relaxed by taking into account a kinetic 
interface undercooling and deviations from chemical equi- 
librium at the alloy's solidification front fl, • 

A number of models 2, 9, 10, 11, 12, 13] have been pro- 
posed to account for solute trapping and related phenom- 
ena observed during rapid phase transformations. One of 
the well-established boundary conditions for solute redis- 
tribution can be taken from the continuous growth model 
(CGM) applied to solute trapping by Aziz and Kaplan 
[3, [l3, [l^. The CGM assumes alloy solidification at a 
"rough interface" ; i.e, all interface sites are potential sites 
for crystallization events. With a high solidification rate, 
the atom can be trapped on a high-energy site of the 
crystal lattice. This leads to a local nonequilibrium on 
the interface and to the formation of metastable solids 
(see examples in Ref. [11] )■ As a result, the solute parti- 
tioning function at the solid-liquid interface is described 
by Wm 



k{V) = 



ke + V/Vl 



DI 



1 + V/Vi 



DI 



(4) 



where Vdi is the speed of diffusion at the interface and 
ke is the value of the equilibrium partition coefficient 
given by Eq. ([T]), i.e., with the negligible interface mo- 
tion, V ^ 0. Equation ^ evaluates the ratio Xs/Xl 
at the interface for dilute solutions of B ( "solute" ) in A 
("solvent"). 

The interfacial diffusion speed Vdi is the kinetic pa- 
rameter describing the deviation from chemical equilib- 
rium at the interface. It has been defined as the ratio 
between the diffusion coefficient Dj at the interface and 
the characteristic distance A for the diffusion jump P.[l3j: 
Vdi = Di/X. The distance A is assumed to be equal to 
the width of the solid-liquid interface (few interatomic 
distances) and the diffusion jumps are taken along the 
direction of growth. Therefore, this definition for Vdi 
is corrected by results of molecular dynamic simulations 
fl6']. They include diffusion in all spatial directions; i.e., 
the diffusion speed is Vdi = 6Dj/X, where the factor 
of 6 accounts for the possibility of jumps along the six 
{±x, y, z) Cartesian axes. 

Outcomes following from the solute partitioning func- 
tion ([4]) were compared in the modeling of solute trap- 
ping using numerical computations based on the phase- 
field theory of alloys solidification. Wheeler et al. [l3] 
naturally included an energy penalty for high composi- 
tion gradients in the liquid that supresses the partition- 
ing of solute at a rapidly moving interface and leads to 
solute trapping. They also showed that the construc- 
tion of common tangents to the curves of free energy (in 
the spirit of Baker and Cahn [l§]) has to be defined for 
nonequilibrium concentrations which already depend on 
the solidification velocity. In order to eliminate or re- 
duce the solute trapping effect by the diffuse interface at 
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small growth velocity, Karma and co-workers proposed 
an ad hoc suitable antitrapping condition to the diffu- 
sion flux 19]. These works [13, [I3| showed that when 
the solute trapping effect comes to modeling alloy so- 
lidification with both phase and concentration fields, a 
crucial issue arises concerning the relative magnitudes of 
the gradients of the two fields within the solidification 
front as well as the relative thickness of the concentra- 
tion jump interface. Additionally, Conti poj investigated 
the usual one-dimensional (ID) formulation of the phase 
field model without the concentration gradient correc- 
tions of Wheeler et al. . He resolved the governing 
equations numerically for the interface temperature and 
the solute concentration field as a function of the growth 
velocity. The partition coefficient k{V) is monotonically 
increasing towards unity at large growth rates following 
the predictions of the continuous growth model (|4]) . How- 
ever, in contrast to the results of natural experiments 
P, [E [3, [E [1] , numeric predictions [l7, 20] were not able 
to reach the complete chemically partitionless (diffusion- 
less) solidification at a finite solidification velocity. 

One of the deficiencies of the function ^ is the diffi- 
culty to describe complete solute trapping at the finite 
solidification velocity: Equation ([4]) predicts k ^ 1 only 
with V oo. Contrary to this prediction, a transition 
to partitionless solidification occurs at a finite solidifica- 
tion velocity as it has been shown in numerous experi- 
ments [H, [1, [3, [S [sl • Molecular dynamic simulations also 
show that the transition to complete solute trapping is 
observed at a finite interface velocity in rapid solidifica- 
tion of a binary system [l6j . Therefore, as an extension 
of Eq. (|4|) , a generalized function for solute partitioning 
in the case of local nonequilibrium solute diffusion within 
the approximation of a dilute system has been introduced 
by Sobolev [Hf. This yields 



HV) 



(1 - vyv^)k, + V/Vi 



DI 



1 - v^/vl + V/Vi 



V <Vi 



D, 



DI 



kiV) = 1, V> Vd. 



(5) 



The diffusion speed Vd introduced in Eq. ^ is the char- 
acteristic bulk speed. It is defined as a maximum speed 
for the solute diffusion propagation or as a speed for the 
front of solute diffusion profile. In particular, the speed 
Vd is obtained by the speed of propagation of the plane 
harmonic wave away from the solid-liquid interface (see 
the Appendix in Ref. [13] )■ As the velocity V of the 
interface is comparable by magnitude to the speed Vd, 
the high frequency limit takes place: ujtd » 1, where 
Lu is the real cyclic frequency of the plane harmonic wave 
and td is the time for relaxation of the diffusion flux to 
its steady state. In this case, Vd has to be considered 
finite and it is defined as Vd = {D/tdY^'^, where D is 
the diffusion coefficient in bulk liquid. 

In the local equilibrium limit, i.e., when the bulk diffu- 
sive speed is infinite, Vd ^ oo, expression © reduces to 
the function k(V) which takes into account the deviation 
from local equilibrium at the interface only as described 



by Eq. ([4|). The function ([5]) includes the deviation from 
local equilibrium at the interface (introducing interfacial 
diffusion speed Vdi) and in the bulk liquid (introducing 
diffusive speed Vd in the bulk liquid). As Eq. (O shows, 
complete solute trapping k{V) — 1 proceeds at = Vd- 
This result has been introduced by Sobolev from a postu- 
lation about the zero value for the diffusion coefficient at 
V > Vd- The next section further details that the con- 
dition for complete solute trapping follows directly from 
the analysis of solute diffusion flux. 



III. DIFFUSION MASS TRANSPORT AND 
SOLUTE TRAPPING 

In ID solidification along the z-axis, the mass balance 
is given by 



dx 



dj_ 



(6) 



where t is the time and J is the diffusion flux. To consider 
solute trapping in ID local nonequilibrium solidification, 
we take one of the results from a model of rapid phase 
transitions [1^ . Using this model, the evolution equation 
for diffusion flux J along the z-axis is described by 



J ^ M 



d_ f ds_ 

dz \dX 



' dz'^ 



dJ_ 

' dt 



(7) 



where s is the entropy density, the factor proportional 
to the correlation length, and M is the diffusion mobility 
of atoms. The latter is defined by 



M^TD, D = D{d{A^i)/dXy 



(8) 



where D is the diffusion coefficient and A/i is the differ- 
ence of chemical potentials between solvent and solute. 
From known thermodynamic expressions [23] one can ac- 
cept that 



ds 
dX 



TD djA^i) 
TD dX 



^, (9) 
TD ^ ' 



where td is the time for diffusion ffux relaxation to its 
steady state. Then, omitting the term responsible for 
atomic correlation, i.e., assuming that Ex = 0, one can 
get from Eq. ([7|) the expression 



J = M 



J_a(A^ _ TD dJ 
T dz 'Wdt 



(10) 



Using Eq. 
follows 



the evolution equation (fTII|) results as 



J 



-D 



^{A^l) 
dz 



TD 



dJ_ 
'dt' 



(11) 



We find the solution for the diffusion flux J which has 
signiflcance in the analysis of solute trapping. Using the 
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expression for the diffusion speed, Vd 
Eqs. (O and PT|) one gets 



9V 



dJ_ 
'dt 



D 



dz'^ 



{D/tdY''^, from 



(12) 



Equation (jl2p is a partial differential equation of hyper- 
bolic type. It describes the flux J in the so-called "hy- 
perbolic evolution" , which proceeds with a sharp front of 
the profile for the solute transport. It occurs due to both 
the diffusive and propagative nature of the transport in 
the high frequency limit ujtd >> 1 with V ^ Vd- 

For a steady-state regime of interfacial motion Eq. (fT^ 
takes the form 



D 



dz 



0, 



(13) 



which is true in a reference frame moving at constant 
velocity V with the interface z = 0. A general solution 
of Eq. 1131) is 



J{z) 



ci + C2 exp 



Vz 



D{1 VyVl,) 



(14) 



To define a particular solution one can assume the fol- 
lowing boundary conditions: the balance on the interface 
is J{z = 0) = V{X*^ - Xg), and the flux is hmited by 
the expression J(z — > oo) = far from the interface with 
z — > oo. The latter condition gives Ci — Q for any veloc- 
ity V and one gets C2 — for V > Vd- Also, from the 
interfacial balance with z = one gets C2 = V{X]^ — Xg) 
for V < Vd- As a result, solution (fT4|) transforms into 
the particular solution 



J{z) 
^V{Xl 
J{z 



X*g) exp 



Vz 



D{1 



v^ivl) 



0, 



V>Vd. 



V <VD: 

(15) 



where X£ and Xg are the liquid concentration and solid 
concentration, respectively, on the interface. 

Solution (fTS]) gives the condition for complete solute 
trapping with finite velocity V > Vd- This is expressed 
by the expression for the solute partitioning function 



k{V) 
k{V) 



XslXl ^ 1, 
1, 



X%^ X 



V<Vd, 
V > Vd- 



(16) 



The latter condition in Eq. defines the equality of 

the concentrations in the phases and leads to complete 
solute trapping. 

To obtain an explicit form for the solute partitioning 
function p6p . we analyze the balance of diffusion fluxes 
on the interface. Taking again the steady state regime 
of solidification constant velocity V, from the system (jS]) 
and (jlip one can obtain the equations 



dJ 

dz 



V- 



.dX 
dz ' 



J 



-D 



diAn) 
dz 



tdV 



dJ 

dz ' 



(17) 



from which we get the single equation for the diffusion 
flux J. This yields 



J 



djAfi) 
dz 



^V^d_X 
Vl dz 



-JT- (18) 



The above-defined thermodynamic parameter D and the 
diffusion speed Vd in bulk have been taken into account. 
Defining the gradient of the difference of the chemical 
potentials as d[A^j)/dz = [d{A^i)/dX]dX/dz, Eq. pS]) 



gives 



1 



d{An) 
dz 



(19) 



This equation is a general expression for the steady dif- 
fusion flux into the liquid from the interface. Within 
the local equilibrium limit Vd oo one can obtain the 
known Fickian approximation which has been used pre- 
viously for analysis of solute trapping [25, 26]. 

Analytical solutions [13] for solidification under local 
nonequilibrium diffusion show that the concentration in 
both phases becomes equal to the initial (nominal) con- 
centration and the diffusion flux is absent for V > Vd- It 
is also given by Eq. (fT5| . Therefore, in addition to Eq. 
(Uni), one can finally obtain 



d(A^ 

dz 



V <V, 



J = 0, 



V > Vd- 



(20) 



At the phase interface one assumes in Eq. (|20p. first, 
that the term D{{dAii) / dX)^^ is proportional to con- 
centration such that 



DiXl 
RT ' 



(21) 



Second, the chemical inhomogeneity (solutal segregation) 
exists due to the jump of the chemical potential A/i 
which has the interfacial gradient —d{AiJ,)/dz = Afi/Wo 
at a small distance Wq of the order of a few inter- 
atomic distances. Third, in an approximation of ideal 
(or even real) solutions, one can assume for the interfa- 
cial difference of chemical potentials A// = A/i^(X) — 
Afi^iX) = RT{lnXl - \nXl) - i?T(lnX| - \nXt;) = 
RT\n{X*g/Xl) - RT\n{Xl/Xl) = RT[\iYk{V) - In fee], 
where X£ and Xg are equilibrium concentrations on the 
interface from the liquid phase and solid phase, respec- 
tively, and k{V) is the ratio of concentrations on the in- 
terface defined by Eq. pB)) . Taking into account these 
last evaluations, one can get for the chemical potential 
gradient the expression 



d{An) _ RT k{V) 
dz Wo ke 



(22) 



Integration of the balance (fT7)) on the interface gives 
the flux 



J = V{Xl - X*g). 



(23) 
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Substituting Eqs. (|2T|) and ((22)) into the expression for 
diffusion flux Eq. ([^0]) with using the balance ([^5]) gives 
the following expression for the solute partitioning func- 
tion: 



1 



In 



k{V) 



V 



DI 



[l-k{V% 



V <Vd 



k{V)=X*s/Xl^l, V>Vd, 



(24) 



where Vdi — Di/Wq is the speed for solute diffusion on 
the interface. 

Equation ((24|) gives the evaluation of solute trapping 
effect through the solute partitioning function k{V) de- 
rived initially from the analysis of the evolution equation 
([7]) for the diffusion flux J. This equation takes into ac- 
count finite diffusion speeds on the interface and in bulk 
liquid. The introduction of these two speeds is a conse- 
quence of the local nonequilibrium both on the interface 
and in bulk liquid. As Eq. (|24|) shows, the complete so- 
lute trapping k{V) — 1 proceeds a.t V — Vd- Equation 
(I24p transforms into a previously known expression for 
the function k{V) derived in Refs. [25, 26] with relaxing 
local equilibrium on the interface and using local equi- 
librium in bulk liquid {Vd — *■ oo) for the diluted binary 
system {XI << 1). 



IV. SOLUTE PARTITIONING FUNCTION 

We use a model of diffusion in which particles move 
by diffusion jumps in random time between two phases 
(states). This model was called the "two- level model of 
diffusion" and it was introduced in the context of various 
application, e.g., in chromatography [28| or for a longi- 
tudinal solute dispersion in a tube with flowing water 
(Taylor's dispersion) [2^. 

Let Pi{t, z) be the probability density of a particle po- 
sition in the phase i — L ov va. the phase i — S aX the 
moment t. Then local conservation of the probability 
density in a point with coordinate z belonging to the 
phase i is defined by 



dt 



dz 



(25) 



If the interface moves with a velocity comparable to the 
solute diffusion speed Vd in bulk phases, then the flux 
Ji{t,z) of the density probability depends on the pre- 
history of the diffusion process. The flux, therefore, is 
defined by 



Mt, z) 



f D,{t~f 

J —OO 



^ dPdt\z) 
dz 



dt* 



(26) 



The relaxation function Di{t~t*) can be chosen in the 
form Di{t — t*) = Di{0) exp[— (t - t*)/TD] of exponential 
decay. In such a case, Eq. (pS)) is reduced to the Maxwell- 
Kattaneo equation 



at 



A(0) 



dz 



0. 



(27) 



It is accepted in Eq. ([27| that Di (0) is the diffusion coeffi- 
cient at the final moment of relaxation prehistory so that 
Di{Q) = D. System ([25]) and ([27]) gives a single equation 
of a hyperbolic type for the density of probability 



TD- 



or for the flux, 



td 



m 

dt 



dt 



D. 



D 



d^ 

dz^ ' 



dz-^ 



(28) 



(29) 



As was shown in Ref. [3Q|], the density of probability 
described by Eq. ([28]) gives a positive entropy production 
for the particle exchange between two levels (between two 
subsystems or phases). 

Integration of Eq. (j28p by an inflnitesimal layer in- 
cluding an interface leads to the balance 



m 

dz 



td- 



d {VR,) 
dt 



0. 



(30) 



In the steady-state regime one can get the following equa- 
tion for the i-th phase 



dR 



td- 



d{VPi) 



dt 



Ji = — tdV -J— 

dz dz 



V' 



m 

dz 



Ji. 



(31) 



which is true in a reference frame moving with constant 
velocity V and placed on the interface where the balance 
([25)1 is described as dJ^/dz = VdP^/dz. Using Eq. |3T|) . 
the balance (|30l) is 



Jl-Js 



Dl{1- 



dPL 

dz 



^DS 



dPs 
dz 



(32) 



In Eq. (l32|) we introduce the speeds Vdl and Vds 
of interfacial solute diffusion from the liquid and solid 
phases, respectively. They are defined by 

Vdl = Dl/Id = vlId, Vds = Ds/Id = i^s^d, (33) 

where Dl and Ds are the diffusion coefficients in the 
phases, Id scales for diffusion within which the diffusion 
jumps occur in phases (or on the interface), and vl and 
vs are the frequencies of diffusion jumps in phases (or 
on the interface). From the theory of the transitive state 
[3l| one can define the frequencies of atomic jumps as 



= vo exp 



us = vq exp 



RT 

Qd + Am' 

RT 



(34) 
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where i/q is the attempt frequency of atomic jumps of the 
order of the vibrational frequency [1, [12] , Q d the acti- 
vation barrier for atomic diffusion through the interface, 
and A^' is the difference of chemical potentials defined 
by Eqs. ([2]) and Q. Obviously, interfacial equilibrium 
exists for 

i^s/tyL = exp[-AM7(i?T)] = fee. (35) 
From the interfacial balance ([32|) there follows 

d 



Jl-Js 



1 



[DLPL-Ds{V)Ps]AiQ) 



where Ds{V) is the function of the interfacial velocity V 
defined by 



Ds{V)^Ds 
0, 



DS 



Vds «Vdl, 



Vds ~ Vol, 



(37) 



Ds{l - V^IVIj}- 



The function ([57]) describes the following cases: (a) 
Vds 0, negligible diffusion in solid {Ps = 0) in com- 
parison with the diffusion in liquid; (b) Vds ~ Vdl, ap- 
proximate equality for diffusion speeds in the liquid and 
solid around the interface; and (c) Vos oo, condition 
of local equilibrium in the diffusion field of the solid (that 
occurs with high frequency jumps of atoms in solid). 

From now on, the above case (b) for approximate 
equality of diffusion speeds in phases around the inter- 
face is taken. First, we use the finite difference —dx = Id 
in the balance ([5^ . Second, we take into account that 
the factor (f — V^/Vj^) is related to the bulk diffusion. 
Finally, using the definition the balance ([5^ is de- 
scribed by 



Jl — Js 



1--^ ) [VdlPl-VdsPs], 



(38) 



where Vd is the solute diffusion speed in bulk liquid 
around the interface. Using Eqs. ([33| - (p5)) . this balance 
can be rewritten as 



Ji - Js = Vol 1 



D 



kePs]- 



(39) 



For the concentrated binary system the probabilities 
Pl and Ps in Eq. (I39p are directly proportional to the 
atomic concentrations in phases. This leads to 

PL^Xs{i-XL)/n, Ps = XL{i~Xs)/n, (40) 



where f2 is the atomic volume. Therefore, Eq. ([39]) can 
be rewritten as 



D 



:[Xsil - Xl) - keXLil - Xs)] 



VPL 

n 



(41) 



We further use the already obtained result (I20p according 
to which the diffusion flux is absent at V > Vd- Then, 
the difference (j4T|) of fluxes on the interface takes the 
form 



Jl ~ Js - [ ^ - y2 



b 



x[X5(l-XL)-fceXL(l-Xs)]^, V<Vd, 



Jl = Js, V > Vd. (42) 
The net flux (|42ll must be equal to the diffusion flux 

V 



Jd - {Xl - X, 



s) 



n' 



(43) 



From the equality of Eqs. (|^^ and one gets 



V / V2 



x[Xs{l-XL)-keXLil-Xs)], V<Vd, 



Xl^Xs, V>Vd, 



(44) 



in which Vd/ = Vdl is the diffusion speed on the inter- 
face from the liquid phase. Equation (|44p can be eas- 
ily resolved regarding the function fc(V) = Xs/Xl of 
nonequilibrium solute partitioning. This yields 

k{V,Xl) 



(f - vyv^)k, + V/V, 



DI 



[1 - vyv^)[i - (1 - ke)xi] + v/Vdi 

k{V,Xl = Xa) = l, V>Vd, 



, V<Vd, 



(45) 



where X^ is the solute concentration in the liquid at the 
interface. 



DISCUSSION AND COMPARISON WITH 
EXPERIMENTAL DATA 



Expression ([45]) gives the general functional depen- 
dence of solute partitioning at the phase interface for 
concentrated binary systems, and with application to 
rapid solidiflcation, it exhibits the two known limits. In 
the first limit, when solidification proceeds with local 
nonequilibrium at the interface only, i.e., with Vd — > oo, 
Eq. (|45p leads to the solute partitioning function of Aziz 
and Kaplan In the second limit, as the concentration 
X£ of the second dissolved component becomes small, 
i.e., the term (1 — fce)X£ might be negligible in compar- 
ison with the unity, Eq. (I45p transforms into Eq. ([5]) as 
suggested by Sobolev ^] . 
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FIG. 1: Predictions of the model given by Eq. (|47)| . Con- 
stants of the binary system are equilibrium partition coeffi- 
cient fee = O.f, bulk diffusion speed Vd ~ 25 (m/s), and in- 
terface diffusion speed Vdi ~ 20 (m/s). The curves present: 
diluted system, (f — fce)Xo << 1 (dotted line) ; slightly concen- 
trated system, Xo = O.fO (dashed line); concentrated system 
Xa = 0.25, (dash-dotted line); and equiconcentrated system, 
Xo = 0.50 (solid line). 



From the analytical solution of the problem of rapid 
solidification under the steady-state regime ^2T\ , the con- 
centration at the planar interface is given by 



k{vy 



X*L — Xq, 



V <Vd, 



V>Vd, 



(46) 



where Xq is the nominal (initial) concentration of the 
solute in the system. In accordance with the solutions 
obtained in Refs. [l^j, a source of concentration pertur- 
bations, i.e., the solid-liquid interface, moving at a ve- 
locity V equal to or higher than the maximum speed Vd 
of these perturbations, cannot change the concentration 
or create the concentration profile ahead of itself. As 
a result for the interface, one obtains in Eq. (|46| that 
X]^ — Xg — Xq with V > V]j. Then the substitution of 
Eq. (j46l) into Eq. (|45| leads to the following expression 
for nonequilibrium solute partitioning function: 



k{V,XQ) 



(1 - VyV^)[ke + (1 - ke)Xo] + V/Vl 



DI 



1 - V^lVl + VjVi 



'DI 



k(y,XQ)^\, v>Vd. 



V <Vd, 



(47) 



Figure[T]demonstrates the behavior of solute partition- 
ing, Eq. (|^7)l . as a function of the interface velocity at 



FIG. 2: Nonequilibrium solute partitioning function k{y,Xo) 
given by the various models. Constants of the binary system 
are nominal concentration of a solute Xq = 0.05 mole frac- 
tion, equilibrium partition coefficient ke — 0.22, bulk diffusion 
speed Vd = 19 (m/s), and interface diffusion speed Vdi = 16 
(m/s). The dotted line is given by the model of Aziz [T^ for 
the diluted system (1 — ke)Xo << 1, the dashed line is given 
by the model of Aziz and Kaplan 01, the dash-dotted line is 
given by the model of Sobolev [Sfj] for diluted system, and the 
solid line is predicted by the present model given by Eq. ()47p . 



various nominal solute concentrations. As the system 
deviates from a diluted one, the trapping of a solute be- 
comes much more pronounced. Also, Eq. (j47p shows that, 
independently from the solute concentration within the 
system, the complete solute trapping k{V,XQ) = 1 pro- 
ceeds when the interface velocity becomes equal to or 
greater than the diffusion speed, i.e., with V >Vd- The 
condition of equality of concentrations in the liquid and 
solid [see Eqs. (|44|) and (|45|) ] means that the lines of the 
nonequilibrium kinetic liquidus and solidus in the kinetic 
phase diagram are merging. It can also be considered as 
the characteristics of diffusionless processes. 

As a general outcome, Eq. (|47| includes the follow- 
ing important cases for nonequilibrium phase transfor- 
mations: (i) the dilute limit described by Aziz's model 
Eq. ©, 

(1 — ke)XQ << 1, and Yd oo, 

(ii) the dilute limit described by Sobolev's solute parti- 
tioning function, Eq. ([5]), 

(1 — fce)Ao << 1, and with the finite Vd, 

(Hi) the concentrated system described by Aziz and Ka- 
plan's model, Ref. 

Vd oo, for arbitrary concentration Xq. 
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FIG. 3: Solute partitioning versus interface velocity for exper- 
imental data [Sa . 113] on solidification of Si- As alloys. Curves 
I' and 2' are given by Eq. ^} with Vd —* oo for 4.5 at.% 
and 9.0 at.% of As in Si, respectively (that gives the model of 
Azis and Kaplan 01). They describe experiment at small and 
moderate solidification velocities. Curves 1 and 2 are given by 
(|47ll with the finite speed Vd for 4.5 at.% and 9.0 at.% of As 
in Si, respectively. These show ability to describe experiment 
in a whole region of investigated solidification velocities for 
both alloys. Data for calculations are given in Table U 



In comparison with the present model's prediction de- 
scribed by Eq. ((47|) these limits are plotted in Fig. [2l 

Figure [3] exhibits theoretical predictions for solute par- 
titioning in comparison with experimental data on the 
solidification of Si-As alloys. Introducing the deviation 
from equilibrium at both the interface and bulk liquid al- 
lows one to describe the whole set of experimental data. 
Particularly, the complete solute trapping is predicted by 
Eq. ^ for Si-4.5 at.%As with Vd = 2.5 m/s and for 
Si-9.0 at.%As with Vd = 2.1 m/s (Tabled. This pro- 
vides a much better agreement with experiments than 
that shown by the Aziz-Kaplan model. 

As can be seen in Fig. [31 predictions of the model of 
Azis and Kaplan [Eq. (|47| with Vd — > oo] disagree with 
experimental data in the region 1.7 < V (m/s)< 2.2 of 
solidification velocities. One may note that at the same 
solidification velocity, i.e., below about V — 2 (m/s), the 
"interface temperature - velocity" relationship also ex- 
hibits a clear deviation from experimental data (see Fig. 
11 in Ref. ^3?!). One may also attribute this deviation 
to the increasing influence of local nonequilibrium solute 
diffusion around the interface and intensive solute trap- 
ping. Thermodynamic analysis and numeric evaluations 
confirm the idea about the pronounced influence of local 
equilibrium in bulk liquid on solute trapping and "inter- 
face temperature - velocity" relationship at high solid- 
ification velocity p^ . [s^. This example confirms that 
local nonequilibrium in the solute diffusion field is re- 
sponsible for nonequilibrium effects appearing in rapid 



solidification (such as solute trapping and solute drag) 
and essential influence on the interface response func- 
tions (temperature, concentration, velocity) [ssj . Thus, 
the agreement between Eq. (|47p and experimental data 
demonstrates the pronounced effect of deviation from lo- 
cal equilibrium in bulk liquid on solute trapping at higher 
solidification velocity. 

Summarizing the behavior for solute partitioning 
shown in Figs. [1][31 one can conclude that during rapid 
solidification the consequences of deviations from local 
chemical equilibrium are threefold. First, the partition 
coefficient becomes dependent on the growth velocity. 
Second, the liquidus and solidus lines approach each 
other. For these two cases it can be enough to introduce 
into the theory deviation from local equilibrium at the 
interface only. Third, in the extreme case (if the solid- 
ification velocity is equal to or greater than the atomic 
diffusive speed in bulk liquid) the partition coefficient 
k{V) becomes unity and the liquidus and solidus lines 
coincide. This leads to a solid being far from chemical 
equilibrium upon diffusionless solidification. Such three 
conditions are of special importance in the preparation 
of metastable supersaturated solutions p^ . 



VI. CONCLUSIONS 

Solute trapping in rapid solidification of a binary al- 
loy's system has been considered. It has been shown that 
the condition for complete solute trapping leading to dif- 
fusionless solidification follows directly from the solution 
for the diffusion task. This task assumes both the low- 
frequency regime (purely diffusion) and high-frequency 
regime (diffusion and propagative regime) of atomic mo- 
tion in a phenomenological statement. 

The two-level model has been used to define the solute 
partitioning function. This model has been used previ- 
ously (e.g., in chromatography and for investigation of 
longitudinal solute dispersion), and it has been formally 
reduced to expressions for an extended version of the con- 
tinuous growth model. The extended version adopts two 
kinetic parameters: solute diffusion speed Vdi on the in- 
terface and solute diffusion speed Vd in bulk liquid. 

A condition of complete solute trapping at the finite so- 
lidification velocity equal to the diffusion speed, V = Vd, 
has been found. This fact is expressed by the general ex- 
pression p6)) for the solute partitioning function. This 
condition defines the equality of the concentration in the 
phases and describes complete solute trapping. Analysis 
leads to concrete forms for the solute partitioning func- 
tion. The first function is given by Eq. (|24p and the 
second function for solute partitioning is described by 
Eq. (|^5)l . Both these functions predict a sharp finishing 
of solute trapping and the onset of diffusionless crystal 
growth at the solidification velocity V equal to the solute 
diffusion speed Vd in bulk liquid. A concrete expression 
for the liquid concentration at the interface allows us 
to give predictions comparable with experimental data. 
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TABLE I: Interface difFusion speed Vdi and bulk diffusion speed Vd for binary systems used in the calculations of the partitioning 
function k{V,Xo) at the solid-liquid planar interface. Equilibrium partition coefficient is taken as ke=0.3 from Refs. [33l. |3^. 





JDilldiy bybLcili 


Vdi l,m/sj 


Vd l.m/sj 




Aziz cHid Kfl.pl&ii's inodBl, R,6f. ^] 




0.37 


_ 


m 


Aziz and Kaplan's model, Ref. [2] 


Si - 9 at.% As 


0.46 
0.37 




[33] 

m 


Sobolev's solute partitioning function, Eq. ([Sjl 


Si - 4.5 at.% As 
and Si - 9 at.% As 


0.75 


2.7 


[21] 


Present model, Eq. (|T7|l 


Si - 4.5 at.% As 


0.8 


2.5 


current data 




Si - 9 at.% As 


0.8 


2.1 


[22] 



The model predicts the complete behavior for the so- 
lute partitioning function dependent on the solidification 
velocity and alloy concentration. In comparison with the 
experimental data of Aziz et al. on solidification of Si- 
As ahoys [M.J. Aziz et al., J. Cryst. Growth 148, 172 
(1995); Acta Mater. 48, 4797 (2000)], the model weU 
predicts deviation of the solute partitioning from equilib- 
rium and complete solute trapping (Fig. [3]). The tran- 
sition from chemically partition growth to diffusionless 
growth at V = Vd occurs sharply. As has been shown 
for dendritic growth [s^ such a sharp transition leads 
to an abrupt exchange of growth kinetics in consistency 



with experimental data. 
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Numerous experimental data on the rapid solidification of binary systems exfiibit tlie formation 
of metastable solid phases with the initial (nominal) chemical composition. This fact is explained 
by complete solute trapping leading to diffusionless (chemically partitionless) solidification at 
a finite growth velocity of crystals. Special attention is paid to developing a model of rapid 
solidification which describes a transition from chemically partitioned to diffusionless growth of 
crystals. Analytical treatments lead to the condition for complete solute trapping which directly 
follows from the analysis of the solute diffusion around the solid-liquid interface and atomic 
attachment and detachment at the interface. The resulting equations for the fiux balance at 
the interface take into account two kinetic parameters: diffusion speed Vdi on the interface and 
diffusion speed Vd in bulk phases. The model describes experimental data on nonequilibrium solute 
partitioning in solidification of Si-As alloys [M.J. Aziz et al., J. Cryst. Growth 148, 172 (1995); 
Acta Mater. 48, 4797 (2000) 1995-2000] for the whole range of solidification velocity investigated. 
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I. INTRODUCTION 

A concept of "solute trapping" has been introduced 
to define the processes of solute redistribution at the in- 
terface which are accompanied by (i) the increasing of 
chemical potential :lll|, and (m) the deviation of the parti- 
tion coefficient for solute distribution towards unity from 
its equihbrium value (independently of the sign of the 
chemical potential) 

In experimental investigations of rapid solidification, a 
complete solute trapping leading to diffusionless (chem- 
ically partitionless) solidification has been first observed 
by Hultgren et al. and Duwez et al. in experiments on 
rapid solidification 0]. They showed that rapidly solid- 
ifying alloy systems lead to the originating of supersat- 
urated solid solution with the initial (nominal) chemical 
composition of the alloy. Later on, crystal microstruc- 
tures with the initial chemical composition were found 
by Chalmers et al. in rapidly solidified pre-dendritic and 
dendritic patterns [1]. 

Backer and Cahn [l| have shown that with the finite 
solidification velocity in a Cd-Zn system the coefficient 
of Cd distribution becomes equal to the unit that char- 
acterizes diffusionless solidification. This fact has been 
confirmed in many binary systems by Miroshnichenko 
. He investigated dendritic crystal microstructure after 
quenching from the liquid state by splat quenching and 
melt spinning methods. Results of Miroshnichenko's mi- 
crostructural analysis show that at a cooling rate greater 
than some critical value (depending on an alloy and ex- 
perimental method this value is in the range 10^ — 10^ 
K/s) a core of main stems of dendrites has initial (nom- 
inal) chemical composition of the alloy. A critical value 
for undercooling in the transition to purely thermally- 
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controlled growth with a homogeneous distribution of 
chemical composition in Ni-B solidifying samples pro- 
cessed by electromagnetic levitation facility has been ob- 
tained by Eckler et al. Finally, it is necessary to note 
that many eutectic systems undergo chemicallypartition- 
less solidification with the initial composition [5] that can 
be explained by the transition to diffusionless solidifica- 
tion 

As a consequence, experimental investigations [l|,[l,0, 
@ show that with increasing driving force of solidifi- 
cation solute traps much more pronounced by solidifying 
microstructure. At the finite value of critical governing 
parameters (undercooling, cooling rate or temperature 
gradient) complete solute trapping occurs. Because the 
finite value of the governing parameter defines the con- 
crete solidification velocity, complete solute trapping and 
diffusionless solidification begin to proceed with the fixed 
critical growth velocity of crystals. 

The main purpose of the present paper is to describe a 
model for solute trapping and transition from chemically 
partitioned to diffusionless solidification in a binary sys- 
tem. Using the local non-equilibrium approach to rapid 
solidification, the analysis of diffusion mass transport in 
bulk phases together with conditions of atomic attach- 
ment/detachment on the solid-liquid interface is given. 

The paper is organized as follows. In Sec. |lTl previous 
investigations on solute trapping are shortly overviewed. 
In Sec. mil the analysis of solute diffusion leading to pro- 
nounced solute trapping and complete solute trapping is 
given. The non-equilibrium solute partitioning function 
for atoms on the interface is derived in Sec. IIVI A com- 
parison with previous models and experimental data on 
solidification of binary systems is presented in Sec. [V] 
Finally, in Sec. IVII conclusions of the work are summa- 
rized. 
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II. PREVIOUS INVESTIGATIONS 



B. Non-equilibrium 



For the simplest case of atomic system, let us consider 
an isobaric and isothermal binary system (the pressure 
P and temperature T are constant) with concentration 
Xa and Xb of atoms A and B, respectively. In this arti- 
cle, we denote X as the concentration of the atoms of B 
sort. For brief overview, we summarize equilibrium and 
non-equilibrium solute distribution on the solid-liquid in- 
terface. 



A. Equilibrium 

In equilibrium, the concentration of atoms X at the 
phase interface is not equal from both sides of the in- 
terface due to different solubility of atoms in phases. 
During the equilibrium co-existence of phases (gas-solid, 
liquid-solid, gas-liquid) the atoms are distributed along 
the interface in consistency with the diagram of a phase 
state. A difference in atomic concentration in phases at 
the interface can be characterized by the equilibrium co- 
efficient fee of atomic distribution between phases. For 
equilibrium co-existence of phases (e.g., between crystal 
and melt, vapor and crystal, crystal and liquid), the co- 
efficient fee can be expressed in the following general form 



ke{XL,Xs,T) 



X 
Xf 



s — 



exp 



(1) 



In Eq. ([T]), A£ and Xg are the mole fractions of the B 
component in the liquid phase (L), or crystal (S), respec- 
tively, R is the gas constant, and Afj,' is the difference in 
chemical potentials described by 



(2) 



with 



^f^B = t^BS - t^BL, ^^J''A^ ^^■'AS~ ^J''ALl (3) 

where A/i^ and A/i'^ are the driving forces for redistri- 
bution of atoms A and B, respectively, which are defined 
by redistribution potentials /i^ and /i'^ for phases L and 
S. The differences A/x^ and A/i^, Eq. ([3]), define the 
sign of A/i' in Eq. ([2]). For instance, if A/i' is negative 
(A/z^ > A/i'g), one has k^. < \ - the case of smaller solu- 
bility of atoms B in the phase S in comparison with their 
solubility in the phase L. 

As a general characteristic of phase equilibria in binary 
systems, expression ([!]), together with Eqs. ^ and ([3]), 
is usually considered as a measure of the driving force 
for atomic redistribution at the phase interface. It can 
also be considered as one of the main parameters for con- 
struction of diagrams of a phase state. 



Expressions ([T])-(l3]) assume local equilibrium at the in- 
terface, which is a useful approximation for many systems 
transforming at small interface velocities. At a large driv- 
ing force for the interface advancing and with increasing 
of the interface velocity, the local equilibrium is not main- 
tained pj]. Therefore, the condition for local interfacial 
equilibrium was relaxed by taking into account a kinetic 
interface undercooling and deviations from chemical equi- 
librium at the alloy's solidification front ^, . 

A number of models [1,1,113,111111111 have been pro- 
posed to account for solute trapping and related phenom- 
ena observed during rapid phase transformations. One of 
the well-established boundary conditions for solute redis- 
tribution can be taken from the continuous growth model 
(CGM) applied to solute trapping by Aziz and Kaplan 
[3, [13, [li|. The CGM assumes alloy solidification at a 
"rough interface" , i.e all interface sites are potential sites 
for crystallization events. With a high solidification rate, 
the atom can be trapped on a high energy site of the 
crystal lattice. This leads to a local non-equilibrium on 
the interface and to the formation of metastable solids 
(see examples in Ref. [Ts']). As a result, the solute parti- 
tioning function at the solid-liquid interface is described 
by ilil 



k{V) = 



ke + V/Vl 



DI 



1 + V/Vi 



DI 



(4) 



where Vdj is the speed of diffusion at the interface, and fee 
is the value of the equilibrium partition coefficient given 
by Eq. i.e. with the negligible interface motion, V — > 
0. Eq. Q evaluates the ratio Xs/X^ at the interface for 
dilute solutions of B ("solute") in A ("solvent"). 

The interfacial diffusion speed, Fp/, is the kinetic pa- 
rameter describing deviation from chemical equilibrium 
at the interface. It has been defined as a ratio between 
the diffusion coefficient, Dj, at the interface and the 
characteristic distance. A, for the diffusion jump [1, [11]: 
Vdi = Dj/X. The distance. A, is assumed to be equal 
to the width of the solid-liquid interface (few interatomic 
distances) and the diffusion jumps are taken along the 
direction of growth. Therefore, this definition for Vdi 
is corrected by results of molecular dynamic simulations 
[l6| . They include diffusion in all spatial directions, i.e., 
the diffusion speed is Vdi — QDi/X, where the factor 
6 accounts for the possibility of jumps along the six 
(±a;, y, z) Cartesian axes. 

Outcomes following from the solute partitioning func- 
tion ([4]) were compared in modeling of solute trapping us- 
ing numerical computations based on the phase-field the- 
ory of alloy's solidification. Wheeler et al. ,17] naturally 
included an energy penalty for high composition gradi- 
ents in the liquid that supresses the partitioning of solute 
at a rapidly moving interface and leads to solute trap- 
ping. They also showed that the construction of common 
tangent to curves of free energy (in a spirit of Baker and 
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Cahn [T^) has to be defined for non-equiUbrium concen- 
trations which aheady depend on soUdification velocity. 
In order to eliminate or reduce the solute trapping effect 
by the diffuse interface at small growth velocity, Karma et 
al. proposed an ad hoc suitable antitrapping condition to 
the diffusion flux [l^ . These works [13, Il9l | showed that 
when the solute trapping effect comes to modeling alloy 
solidification with both phase and concentration fields, 
a crucial issue arises concerning the relative magnitudes 
of the gradients of the two fields within the solidification 
front as well as the relative thickness of the concentration 
jump interface. Additionally, Conti ^^oj investigated ID 
usual formulation of the phase field model without con- 
centration gradient corrections of Wheeler et al. ■ He 
resolved the governing equations numerically for the in- 
terface temperature and the solute concentration field as 
a function of the growth velocity. The partition coeffi- 
cient k(V) is monotonically increasing towards unity at 
large growth rates following the predictions of the con- 
tinuous growth model Q. Howevei\ in contrast to the 
results of natural experiments P, H, Ij, H, @| , the numeric 
predictions [13, 1131 were not able to reach the complete 
chemically partitionless (diffusionless) solidification at a 
finite solidification velocity. 

One of the deficiencies of the function (j4]) is the diffi- 
culty to describe complete solute trapping at the finite 
solidification velocity: Eq. ^ predicts fc — > 1 only with 
V oo. Contrary to this prediction, a transition to 
partitionless solidification occurs at a finite solidification 
velocity as it has been shown in numerous experiments 
P, d, d, d, [1] . Molecular dynamic simulations also show 
that the transition to complete solute trapping is ob- 
served at a finite interface velocity in rapid solidification 
of a binary system [l^]. Therefore, as an extension of 
Eq. ([3]), a generalized function for solute partitioning in 
the case of local non-equilibrium solute diffusion within 
the approximation of a dilute system has been introduced 
by Sobolev [Hj. This yields 



HV) = 



(1 - VyVl)ke + V/Vl 



DI 



V/Vi 



V <Vi 



D, 



DI 



k{V) = 1, V> Vd- 



(5) 



The diffusion speed, Vd, introduced into Eq. ([5]) is the 
characteristic bulk speed. It is defined as a maximum 
speed for solute diffusion propagation or as a speed for 
the front of solute diffusion profile. In particular, the 
speed, Vd, is obtained by the speed of propagation of the 
plane harmonic wave away from the solid-liquid interface 
(see Appendix in Ref. ^2^]). As the velocity, V, of the 
interface is comparable by magnitude to the speed, Vd, 
the high frequency limit takes place: ujtd » 1, where 
Lu is the real cyclic frequency of the plane harmonic wave 
and td is the time for relaxation of the diffusion flux to 
its steady state. In this case, Vd has to be considered 
finite and it is defined as Vd = {D/tdY^'^, where D is 
the diffusion coefficient in bulk liquid. 

In the local equilibrium limit, i.e. when the bulk diffu- 



sive speed is infinite, Vd oo, expression ([5]) reduces to 
the function k(V) that takes into account the deviation 
from local equilibrium at the interface only as described 
by Eq. The function ([5]) includes the deviation from 
local equilibrium at the interface (introducing interfacial 
diffusion speed Vdi) and in the bulk liquid (introduc- 
ing diffusive speed Vd in the bulk liquid). As Eq. ([5]) 
shows, the complete solute trapping, k{V) — 1, proceeds 
at V — Vd- This result has been introduced by Sobolev 
from a postulation about the zero value for the diffusion 
coefficient a.t V > Vd- The next section further details 
that the condition for complete solute trapping follows 
directly from the analysis of solute diffusion fiux. 



III. DIFFUSION MASS TRANSPORT AND 
SOLUTE TRAPPING 

In ID solidification along the z-axis, the mass balance 
is given by 



dX 
'dt 



dJ_ 

dz ' 



(6) 



where t is the time and J is the diffusion flux. To consider 
solute trapping in ID local non-equilibrium solidification, 
we take one of the results from a model of rapid phase 
transitions [23^. Using this model, the evolution equation 
for diffusion ffux J along the z-axis is described by 



d_ 

dz 



ds 
dX 



dJ_ 

dt 



(7) 



where s is the entropy density, the factor proportional 
to the correlation length, and M is the diffusion mobility 
of atoms. The latter is defined by 



M = TD, 



D = D{d{A^i)/dXy 



(8) 



where D is the diffusion coefficient and A/i is the differ- 
ence of chemical potentials between solvent and solute. 
From known thermodynamic expressions [2^ one can ac- 
cept that 



ds 
dX 



TD djA^i) 
TD dX 



(9) 

TD 



where td is the time for diffusion ffux relaxation to its 
steady state. Then, omitting the term responsible for 
atomic correlation, i.e. assuming that Sx = 0, one can 
get from Eq. ([7|) the following expression 



J ^ M 



J_a(A^ _ TD dJ 

T dz ^ TD dt 



(10) 



Using Eq. ([8]), the evolution equation (fTO]) results as 
follows 



j^diA^^ dJ 



(11) 
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We find the solution for the diffusion flux J which has 
significance in analysis of solute trapping. Using the ex- 
pression for the diffusion speed, Vd = (D/td)^^^ , from 
equations © and (fTT|) one gets 



92 J 



TD 



9*2 dt 



dJ _ J 



dz 



(12) 



Eq. (jT2|) is a partial differential equation of a hyperbolic 
type. It describes the flux J in the so-called "hyperbolic 
evolution" , which proceeds with the sharp front of profile 
for the solute transport. It occurs due to both diffusive 
and propagative nature of the transport in the high fre- 
quency limit lotd >> 1 with V ~ Vd- 

For a steady-state regime of interfacial motion Eq. p2|) 
takes the following form 



V 



^dJ 
dz 



= 0, 



(13) 



which is true in a reference frame moving at constant 
velocity V with the interface z = 0. A general solution 
of Eq. (dSl) is 



J{z) = ci + C2 exp 



Vz 



D{i-vyvE) 



(14) 



To define a particular solution one can assume the fol- 
lowing boundary conditions: balance on the interface is 
J(z = 0) = V{Xl - Xt;), and the flux is limited by the 
expression J(z —>■ 00) = far from the interface with 
z ^ 00. The latter condition gives Ci = for any veloc- 
ity V and one gets C2 = for y > Vd- Also, from the 
interfacial balance with z = one gets C2 = V{X'[ — Xg) 



for V < Vd- As a result, solution 
the following particular solution 



J{z) 

= V{Xl-X*s)e^p(^- 
J(z) = 0, V>Vd, 



Vz 



D{i-vyv^) 



IM transforms into 



V<Vd, 
(15) 



where X^ and Xg are the liquid concentration and solid 
concentration, respectively, on the interface. 

Solution (jlSp gives condition for complete solute trap- 
ping with the finite velocity V > Vd- This is expressed 
by the expression for the solute partitioning function 



k{V) 
k{V) 



X*s/Xl ^ 1, 



1. 



xi^ X 



V<Vd, 
V > Vd- 



(16) 



The latter condition in Eq. pB]) defines equality of con- 
centration in the phases and leads to complete solute 
trapping. 

To obtain the explicit form for the solute partitioning 
function p6p . we analyze the balance of diffusion fluxes 
on the interface. Taking again the steady state regime 
of solidiflcation constant velocity, V, from the system © 
and (jlip one can obtain the following equations 



dJ 

dz 



V- 



.dX 
dz ' 



J 



-D 



diAn) 
dz 



tdV 



dJ 

dz ' 



(17) 



from which we get the single equation for the diffusion 
flux J. This yields 



J 



d(A/i) 

dz 



D 



^2 dX 

v^U' 



(18) 



It has been taken into account the defined above thermo- 
dynamic parameter D and the diffusion speed Vd in bulk. 
Defining the gradient of difference of chemical potentials 
as d{An)/dz = {d{Afi)/dX)dX/dz, Eq. HID gives 



V"2 



d(A/x) 
dz 



(19) 



This equation is a general expression for the steady dif- 
fusion flux into the liquid from the interface. Within the 
local equilibrium limit Vd — > 00 one can obtain known 
Fickian approximation that has been used previously for 
analysis of solute trapping [2^ . 

Analytical solutions 27| for solidification under local 
non-equilibrium diffusion show that the concentration in 
both phases becomes equal to the initial (nominal) con- 
centration and the diffusion flux is absent for V > Vd- It 
is also given by Eq. Therefore, in addition to Eq. 

p9)) . one can flnally obtain: 



J = 0, 



V > Vd- 



1 



1^2 



d{Ay) 

dz ' 



V <Vd\ 



(20) 



At the phase interface one assumes in Eq. (|20p. flrst, 
that the term D{{dAfi)/dX)~^ is proportional to con- 
centration such that 



RT 



(21) 



Second, the chemical inhomogeneity (solutal segrega- 
tion) exists due to the jump of the chemical potential 
A/i which has the interfacial gradient —d{A^)/dz = 
A/i/Wo a small distance Wq of the order of few in- 
teratomic distances. Third, in approximation of ideal 
(or even real) solutions, one can assume for the interfa- 
cial difference of chemical potentials: A/j, — Aii^{X) — 
A^^(X) ^ RT{\nXl - lnX£) - RT{\iiX% - InAT^) = 
RT\n{X*s/Xl) - RThi{XI,/Xl) = RT{\nk{V) - In fee), 
where Ar£ and Xg are equilibrium concentrations on the 
interface from the liquid phase and solid phase, respec- 
tively, and k[V) is the ratio of concentrations on the in- 
terface deflned by Eq. pB]) . Taking into account these 
last evaluations, one can get for the chemical potential 
gradient the following expression: 



d{An) _ RT k{V) 



dz 



Wo 



(22) 



Integration of the balance p7|) on the interface gives 
the flux 



J = V{Xl - X*s). 



(23) 
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Substituting Eqs. ((2T|l and ([22)1 into expression for diffu- 
sion flux (pni) with using balance gives the following 
expression for the solute partitioning function 



1 



In 



fee Vz5/ ^ 



k{V) ^ X*s/Xl 



1, 



(24) 



where Vdi = Di/Wq is the speed for solute diffusion on 
the interface. 

Eq. ([M)) gives the evaluation of solute trapping effect 
through the solute partitioning function k(V) derived ini- 
tially from the analysis of the evolution equation ([7]) for 
the diffusion flux J. This equation takes into account 
finite diffusion speeds on the interface and in bulk liquid. 
Introduction these two speeds is a consequence of the 
local non-equilibrium both on the interface and in bulk 
liquid. As Eq. p4)) shows, the complete solute trapping, 
k{V) = 1, proceeds at V = Vd- Eq. ([M)) transforms 
into previously known expression for the function k{V) 
derived in Refs. [1^, |2^ with relaxing local equilibrium 
on the interface and using local equilibrium in bulk liquid 
{Vd oo) for the diluted binary system {X^ << 1). 



IV. SOLUTE PARTITIONING FUNCTION 

We use a model of diffusion in which particles move 
by diffusion jumps in random time between two phases 
(states). This model was called as "two level's model of 
diffusion" and it was introduced in a context of various 
application, e.g., in chromatography [2^ or for a longi- 
tudinal solute dispersion in a tube with flowing water 
(Taylor's dispersion) [2^. 

Let Pi{t, z) is the probability density of a particle po- 
sition in the phase i = L or in the phase i = S' at the 
moment t. Then local conservation of the probability 
density in a point with coordinate z belonging to the 
phase i is defined by 



dt 



dz 



(25) 



If the interface moves with the velocity comparable with 
the solute diffusion speed Vd in bulk phases then the 
flux Ji{t, z) of density probability depends on prehistory 
of the diffusion process. The flux, therefore, is defined by 



D,{t-t*) 



dz 



dt* 



(26) 



Relaxation function Di (t 
D,[t-t*) = A(0)exp(-(t 
cay. In such case, Eq. 
Kattaneo equation 



- 1*) can be chosen in a form 
- t*)/T£i) of exponential de- 
is reduced to the Maxwell- 



at 



A(0) 



dz 



(27) 



It is accepted in Eq. ((27| that -Di(O) is the diffusion coeffi- 
cient at the final moment of relaxation prehistory so that 
Di{Q) = D. System ([25]) and ([27]) gives a single equation 
of a hyperbolic type for the density of probability: 



or for the flux 



TD 



dt 



dt 



d'P. 
dz^ 



= A 



dz'^ 



(28) 



(29) 



As it was shown in Ref. j3]Q], the density of probability 
described by Eq. (^5)1 gives a positive entropy production 
for the particle exchange between two levels (between two 
subsystems or phases). 

Integration of Eq. (pS)) by infinitesimal layer including 
interface leads to the balance 



D.. 



dz 



TD- 



diVI^ 
dt 



= 0. 



(30) 



In the steady-state regime one can get the following equa- 
tion for the i-th phase 



D, 



dP 



td- 



d jVPi, 
dt 



J,. = D, 



dz 



-tdV 



dz 



m 

dz 



Ji 



(31) 



which is true in a reference frame moving with the con- 
stant velocity V and placed on the interface where the 
balance (^5)) is described as dJi/dz = VdPi/dz. Using 
Eq. dSI]), the balance ^ is 



Jl-Js 



Dr. 1 



^DL 



dPL 

dz 



^DS 



dPs 
dz 



(32) 



In Eq. (|32|l we introduce the speeds Vdl and Vds 
of interfacial solute diffusion from the liquid and solid 
phases, respectively. They are defined by 



Vdl = Dl/Id = vlIl 



Vds = Ds/Id = ^sId, (33) 



where Dl and Ds arc the diffusion coefficients in phases. 
Id scale for diffusion within which the diffusion jumps 
occur in phases (or on the interface), and and vs are 
frequencies of diffusion hopes in phases (or on the inter- 
face). From the theory of the transitive state [3l| one 
can define frequencies of atomic jumps as 



= vo exp 



us = vq exp 



RT 

Qd + Am' 

RT 



(34) 
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where vq is the attempt frequency of atomic hopes of the 
order of the vibrational frequency [1, [12] , Q d the acti- 
vation barrier for atomic diffusion through the interface, 
and A/i' is the difference of chemical potentials defined 
by Eqs. ([2]) and Q. Obviously, interfacial equilibrium 
exists for 

i^s/i^L = exp[-A///(i?T)] = ke. (35) 
From the interfacial balance (l32l) there follows 



Jl^Js = -(i- -j^) ^ [DlPl - DsiV)Ps] , (36) 

where Ds{V) is the function of the interfacial velocity V 
defined by 



DsiV) = A 



0, 



Vds «Vdl, 



Vds ~ Vol, 



(37) 



The function ([57)1 describes the following cases: 

(a) Vds — ^ 0: negligible diffusion in solid {Ds = 0) in 
comparison with the diffusion in liquid, 

(b) Vds ~ Vdl' approximate equality for diffusion 
speeds in the liquid and solid around the interface, 

(c) Vds ~^ oo: condition of local equilibrium in the dif- 
fusion field of the solid (that occurs with high frequency 
hopes of atoms in solid). 

From now on, the above case (b) for approximate 
equality of diffusion speeds in phases around the inter- 
face is taken. First, we use the finite difference —dx = Id 
in the balance ([5^ . Second, we take into account that 
the factor (1 — V^^/Vq^) is related to the bulk diffusion. 
Finally, using the definition ((33)) . the balance ((32)) is de- 
scribed by 

Jl-Js=[i- [VdlPl - VdsPs], (38) 

where Vd is the solute diffusion speed in bulk liquid 
around the interface. Using Eqs. ([55)l . ([M)) and (|55)l . 
this balance can be rewritten as 

Jl-Js^ Vdl {^-^^[Pl- kePs] ■ (39) 

For the concentrated binary system the probabilities 
Pl and Ps in Eq. ()39p are directly proportional to the 
atomic concentrations in phases. This leads to 

PL = Xs{i-XL)/n, Ps = XL{i-Xs)/n, (40) 

where D, is the atomic volume. Therefore, Eq. p9p can 
be rewritten as 

/ y2 



D 



:[Xsil - Xl) - keXLil - Xs)] 



Vdl 

n 



(41) 



We further use the already obtained result (I20p according 
to which the diffusion flux is absent at V > Vd- Then, 
the difference (|4T)) of fluxes on the interface takes the 
following form 



T/2 

Jl - JS ^ [ ^ - y2 



D 



,Vdl 



x[Xsil-XL)~keXL{l-Xs)]^, V<Vd, 

Jl = Js, V> Vd. (42) 
The net flux ((42| must be equal to the diffusion flux 

J^^{Xl-Xs)^. (43) 
From equality of Eqs. (|42p and (|^5)l one gets 



V f V^ 



:[Xsil-~XL)-keXLil-Xs)], V<Vd, 



Xl=Xs, V>Vd, 



(44) 



in which Vdi = Vdl is the diffusion speed on the inter- 
face from the liquid phase. Equation (|^^ can be easily 
resolved regarding the function k{V) — Xg/Xj^ of non- 
equilibrium solute partitioning. This yields 



kiV,Xl) 



(1 - VyV^)K + V/Vi 



DI 



(1 - vyv^)[i - (1 - k,)xi] + v/Vdi 

k{V,Xl^Xo)^h V>Vd, 



, V<Vd, 



(45) 



where X^ is the solute concentration in the liquid at the 
interface. 



DISCUSSION AND COMPARISON WITH 
EXPERIMENTAL DATA 



Expression (|45)) gives the general functional depen- 
dence of solute partitioning at the phase interface for 
concentrated binary systems, and, with application to 
rapid solidification, it exhibits the two known limits. In 
the first limit, when solidification proceeds with the local 
non-equilibrium at the interface only, i.e. with Vd ~* oo, 
Eq. (|45p leads to the solute partitioning function of Aziz 
and Kaplan In the second limit, as the concentration 
X£ of the second dissolved component becomes small, 
i.e., the term (1 — fce)X£ might be negligible in com- 
parison with the unity. Eg. ()45p transforms into Eq. ([S]) 
suggested by Sobolev [21|. 
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FIG. 1: Predictions of the model given by Eq. (|47|l . Con- 
stants of a binary system are: equilibrium partition coefficient 
is fee = 0.1, bulk diffusion speed is Vd = 25 (m/s), and inter- 
face diffusion speed is Vdi = 20 (m/s). The curves present: 
diluted system, (1 — fce)Xo << 1, (dotted); slightly concen- 
trated system, Xq = 0.10, (dashed); concentrated system, 
Xo = 0.25, (dashed-dotted); and equi-concentrated system, 
Xo = 0.50, (solid). 



From the analytical solution of the problem of rapid 
solidification under steady-state regime [27] , the concen- 
tration at the planar interface is given by 



X*r 



k{vy 



X*L — Xq-> 



V <Vd, 



V>Vd, 



(46) 



where Xq is the nominal (initial) concentration of the 
solute in the system. In accordance with the solutions 
obtained in Refs. [l^j, a source of concentration pertur- 
bations, i.e. the solid-liquid interface, moving at the ve- 
locity V equals to or higher than the maximum speed Vd 
of these perturbations, cannot change the concentration 
nor create the concentration profile ahead of itself. As 
a result for the interface, one obtains in Eq. (|46| that 
X]^ — Xg = Xq with V > Vjj- Then the substitution of 
Eq. (j46| into Eq. (|45| leads to the following expression 
for non-equilibrium solute partitioning function: 



k{V,Xo) 



(1 - VyV^)[ke + (1 - K)Xo] + V/Vl 



DI 



1 - V^lVl + ViVi 



'DI 



fc(V^,Xo) = l, V>Vd. 



V <Vd, 



(47) 



Figure[T]demonstrates the behavior of solute partition- 
ing, Eq. (|^7)l . as a function of the interface velocity at 



FIG. 2: Non-equilibrium solute partitioning function fe(V, Xo) 
given by the various models. Constants of a binary system 
are: nominal concentration of a solute is Xo — 0.05 mole 
fraction, equilibrium partition coefficient is ke — 0.22, bulk 
diffusion speed is Vd = 19 (m/s), and interface diffusion speed 
is Vdi = 16 (m/s). The dotted line is given by the model of 
Aziz [l^ for the diluted system, (1 — ke)Xo << 1; the dashed 
line is given by the model of Aziz and Kaplan H; the dashed- 
dotted line is given by the model of Sobolev 21] for diluted 
system; and the solid line is predicted by the present model 
given by Eq. (|47t . 



a various nominal solute concentration. As the system 
deviates from the diluted one, the trapping of a solute 
becomes much more pronounced. Also, Eq. (|47| shows 
that, independently from the solute concentration within 
the system, the complete solute trapping, fc(y,Xo) = 1, 
proceeds when the interface velocity becomes equal to or 
greater than the diffusion speed, i.e. with V > Vd- The 
condition of equality of concentrations in the liquid and 
solid [see Eqs. (|44l) and (|45p ] means that the lines of the 
non-equilibrium kinetic liquidus and solidus in the kinetic 
phase diagram are merging. It can also be considered as 
the characteristics of diffusionless process. 

As a general outcome, Eq. (|47)) includes the follow- 
ing important cases for non-equilibrium phase transfor- 
mations: (i) dilute limit described by Aziz's model p^ . 
Eq. Q, 

(1 — ke)XQ << 1, and Vd — > oo, 

(ii) dilute limit described by the Sobolev's solute parti- 
tioning function, Eq. ([5|), 

(1 — ke)XQ << 1, and with the finite Vd, 

(Hi) the concentrated system described by the Aziz and 
Kaplan's model, Ref. 2], 

Vd oo, for arbitrary concentration Xq. 
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FIG. 3: Solute partitioning versus interface velocity for exper- 
imental data [Sa . 113] on solidification of Si- As alloys. Curves 
I' and 2' are given by Eq. ^} with Vd —* oo for 4.5 at.% 
and 9.0 at.% of As in Si, respectively (that gives the model of 
Azis and Kaplan 01). They describe experiment at small and 
moderate solidification velocities. Curves 1 and 2 are given by 
(|47ll with the finite speed Vd for 4.5 at.% and 9.0 at.% of As 
in Si, respectively. These show ability to describe experiment 
in a whole region of investigated solidification velocities for 
both alloys. Data for calculations are given in Table U 



In comparison with the present model's prediction de- 
scribed by Eq. ([47|) these limits are plotted in Fig. [2l 

Figure [3] exhibits theoretical predictions for solute par- 
titioning in comparison with experimental data on solid- 
ification of Si-As alloys. Introducing the deviation from 
equilibrium both at the interface and bulk liquid allows 
one to describe the whole set of experimental data. Par- 
ticularly, the complete solute trapping is predicted by 
Eq. (gTl) for Si-4.5 at.%As with Vd = 2.5 m/s and for 
Si-9.0 at.%As with Vd = 2.1 m/s (Tabled. This pro- 
vides a much better agreement with experiments than 
that shown by Aziz and Kaplan model. 

As it can be seen in Fig. [31 predictions of the model of 
Azis and Kaplan [Eq. (|47| with Vd — > oo] disagree with 
experimental data in the region 1.7 < V (m/s)< 2.2 of 
solidification velocities. One may note that at the same 
solidification velocity, i.e. below about V = 2 (m/s), the 
"interface temperature - velocity" relationship also ex- 
hibits a clear deviation from experimental data (see Fig. 
11 in Ref. [s^l)- One may also attribute this deviation 
to the increasing influence of local non-equilibrium so- 
lute diffusion around the interface and intensive solute 
trapping. Thermodynamic analysis and numeric evalua- 
tions confirm the idea about pronounced influence of local 
equilibrium in bulk liquid on solute trapping and "inter- 
face temperature - velocity" relationship at high solidifl- 
cation velocity [S^, [sHi ■ This example confirms that local 
non-equilibrium in solute diffusion field is responsible for 
non-equilibrium effects appearing in rapid solidification 



(such as solute trapping and solute drag) and essential 
influence on the interface response functions (tempera- 
ture, concentration, velocity) [ssj . Thus, the agreement 
between Eq. (|47p and experimental data demonstrates 
the pronounced elfect of deviation from local equilibrium 
in bulk liquid on solute trapping at higher solidification 
velocity. 

Summarizing the behavior for solute partitioning 
shown in Figs. [1][31 one can conclude that during rapid 
solidification the consequences of deviations from local 
chemical equilibrium are threefold. First, the partition 
coefficient becomes dependent on the growth velocity. 
Second, the liquidus and solidus lines approach each 
other. For these two cases can be enough to introduce 
into the theory deviation from local equilibrium at the 
interface only. Third, in the extreme case (if the solid- 
ification velocity is equal to or greater than the atomic 
diffusive speed in bulk liquid) the partition coefficient 
k{V) becomes unity and the liquidus and solidus lines 
coincide. This leads to a solid being far from chemical 
equilibrium upon diffusionless solidification. Such three 
conditions are of special importance in the preparation 
of metastable supersaturated solutions p^ . 



VI. CONCLUSIONS 

Solute trapping in rapid solidification of a binary al- 
loy's system has been considered. It has been shown 
that the condition for the complete solute trapping lead- 
ing to the diffusionless solidification follows directly from 
the solution for the diffusion task. This task assumes 
both low-frequency regime (purely diffusion) and high- 
frequency regime (diffusion and propagative regime) of 
atomic motion in phenomenological statement. 

The two level model has been used to define the solute 
partitioning function. This model has been used previ- 
ously (e.g., in chromatography and for investigation of 
longitudinal solute dispersion) and it has been formally 
reduced to expressions for extended version of the con- 
tinuous growth model. The extended version adopts two 
kinetic parameters: solute diffusion speed Vdi on the in- 
terface and solute diffusion speed Vd in bulk liquid. 

A condition of complete solute trapping at the finite so- 
lidification velocity equal to the diffusion speed, V = Vd, 
has been found. This fact is expressed by the general ex- 
pression (fT6)) for solute partitioning function. This condi- 
tion defines equality of concentration in the phases and 
describes the complete solute trapping. Analysis leads 
to concrete forms for the solute partitioning function. 
The first function is given by Eq. p4l) and the second 
function for solute partitioning is described by Eq. (|45p . 
These both functions predict the sharp finishing of solute 
trapping and the onset of diffusionless crystal growth at 
the solidification velocity V equal to the solute diffusion 
speed Vd in bulk liquid. The concrete expression for 
liquid concentration at the interface allows to give 
predictions comparable with experimental data. 
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TABLE I: Interface difFusion speed Vdi and bulk diffusion speed Vd for binary systems used in the calculations of the partitioning 
function k{V,Xo) at the solid-liquid planar interface. Equilibrium partition coefficient is taken as ke=0.3 from Refs. [33l. |3^. 



A/Tr\rl*:>l 


jDiiidiy oybLciii 


Vdi l,m/sj 


Vd l.m/sj 




^T~7 o n o taI q "Ki ' c tt^ /^/^ £i^'\ "R iiif rOl 


Ol - ^.O cLL./O -Tib 


n dp, 
0.37 


_ 


Ref. [34] 


Aziz and Kaplan's model, Ref. [2] 


Si - 9 at.% As 


0.46 
0.37 




Ref. [33] 
Ref. [34] 


Sobolev's solute partitioning function, Eq. ([Sjl 


Si - 4.5 at.% As 
and Si - 9 at.% As 


0.75 


2.7 


Ref. [21] 


Present model, Eq. (|T7|l 


Si - 4.5 at.% As 


0.8 


2.5 


current data 




Si - 9 at.% As 


0.8 


2.1 


Ref. [22] 



The model predicts the complete behavior for the so- 
lute partitioning function dependent on solidification ve- 
locity and alloy's concentration. In comparison with ex- 
perimental data of Aziz et al. on solidification of Si- As 
alloys, the model well predicts deviation of the solute 
partitioning from equilibrium and the complete solute 
trapping (Fig. [J]). The transition from chemically par- 
tition growth to diffusionless growth at V — Vd occurs 
sharply. As it has been shown for dendritic growth (36j 
such sharp transition leads to the abrupt exchange of 
growth kinetics in consistency with experimental data. 
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